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Abstract 

We  study  the  problem  of  source  localization  as  a  multiple  hypothesis  testing,  where  each  hypothesis  corresponds  to  the 
event  that  the  source  belongs  to  a  particular  region.  We  use  sequential  hypothesis  tests  based  on  posterior  computations 
to  solve  for  the  correct  hypothesis.  Measurements  corrupted  with  noise  are  used  to  calculate  conditional  posteriors. 
We  prove  that  the  regional  localization  problem  has  geometric  properties  that  allow  correct  detection  almost  surely  in 
the  limit  of  infinite  measurements.  We  present  the  Sense,  Transmit  &  Test  algorithm  that  allows  sequential  sensing, 
communication  and  testing  and  we  analyze  the  accuracy  of  this  distributed  algorithm  and  show  that  the  test  ends  in  a 
finite  time.  We  also  present  numerical  results  illustrating  properties  of  the  suggested  algorithm. 
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1.  Introduction 

1.1.  Problem  description  and  motivation 

Applications  where  source  localization  is  of  great  con¬ 
cern,  vary  between  finding  the  source  of  oil  spills  in  the 
ocean,  determining  cellular  locations,  detecting  an  earth¬ 
quake’s  epicenter,  locating  an  acoustic  source,  or  simply 
finding  an  intruder  in  a  protected  environment.  For  most 
of  these  applications,  it  is  sufficient  to  find  a  region  that 
contains  the  source  rather  than  pinpointing  the  exact  source 
position,  which  relies  most  of  the  time  on  approximations. 

In  this  work  we  consider  the  following  problem:  A 
source  at  an  unknown  location  in  a  bounded  region  Q 
transmits  a  power  signal.  N  sensors  receive  noisy  and 
decayed  versions  of  the  signal,  they  can  communicate  and 
exchange  measurements.  The  environment  Q  is  divided 
into  M  regions  Wa ,  where  a  €  {1, . . . ,  M}.  The  objective 
of  the  sensors  is  to  find  which  region  contains  the  source. 

We  pose  the  problem  as  a  multiple  hypothesis  testing 
problem,  where  hypothesis  Ha  is  true  if  the  source  lies  in 
the  region  Wa.  We  assume  no  prior  knowledge  about  the 
location  of  the  source  and  therefore  model  the  source  loca¬ 
tion  as  a  uniformly  distributed  random  variable  over  the 
environment  Q ,  any  prior  information  about  the  source  lo¬ 
cation  can  be  incorporated  in  the  location  density  function. 
We  adopt  the  log-normal  fading  model  for  the  propagation 
of  the  received  signal  power.  The  noise  added  to  the  log 
of  the  power  is  Gaussian  with  zero  mean  and  a  known 
variance  er2. 
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1.2.  Literature  review 

In  the  classical  source  localization  problem,  a  number 
of  sensors  collaborate  to  locate  the  exact  position  of  a 
source.  The  relation  between  the  position  of  a  source  and 
the  received  signal  strength  (RSS)  is  described  in  (Rap¬ 
poport,  1996;  Proakis  and  Salehi,  2001;  Sayed  et  al.,  2005; 
Chen  et  al.,  2002).  Several  authors  treat  localization  as 
a  nonconvex  optimization  problem  (Hero  III  and  Blatt, 
2005;  Rabbat  and  Nowak,  2004a).  Gradient  descent  al¬ 
gorithms  and  weighted  least  squares  approximations  can 
be  used  to  solve  the  maximum  likelihood  estimation  prob¬ 
lems  but  such  algorithms  tend  to  get  stuck  at  local  opti- 
mas  (Rabbat  and  Nowak,  2004b;  Mao  et  al.,  2007).  Meng 
et  al.  (2008)  approximate  the  nonlinear  nonconvex  opti¬ 
mization  problem  by  a  linear  and  convex  problem.  Hero 
III  and  Blatt  (2005)  use  a  method  of  projection  onto  con¬ 
vex  sets.  A  necessary  and  sufficient  condition  for  the  con¬ 
vergence  of  this  algorithm  is  that  the  source  lies  inside  the 
convex  hull  of  the  sensors.  Properly  placing  the  sensors 
assumes  knowledge  of  the  position  of  the  source. 

Designing  distributed  algorithms  is  in  general  a  prob¬ 
lem  specific  task,  and  many  researchers  from  various  com¬ 
munities  have  looked  at  this  problem.  We  refer  the  reader 
to  (Nedic  and  Ozdaglar,  2009;  Lynch,  1997;  Boyd  et  al., 
2006)  and  references  therein  for  more  details  about  this 
topic. 

The  multiple  hypothesis  problems  are  considerably  more 
difficult  than  the  binary  problem  and  optimality  of  the  pro¬ 
posed  algorithms  is  usually  hard  to  prove.  Some  tests  that 
have  some  asymptotic  optimality  properties  were  devel¬ 
oped  in  the  literature,  but  these  tests  tend  to  be  very  com¬ 
plex  (Savin,  1984;  Baum  and  Veeravalli,  1994;  Armitage, 
1950).  Alternatively  ad  hoc  tests  based  on  repeated  pair¬ 
wise  applications  of  optimal  sequential  hypothesis  tests  (Wald, 
1945)  were  developed  but  these  tests  have  little  optimality 
results,  e.g.,  see  Eisenberg  (1991). 
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1.3.  Contributions 

The  contributions  of  this  paper  are  three- folds. 

First,  we  formulate  the  source  localization  problem  in  a 
novel  multi-hypothesis  testing  setting.  We  analyze  proper¬ 
ties  of  the  Maximum  A  Posteriori  (MAP)  algorithm  that 
requires  the  computation  of  a  finite  number  of  integrals 
which  is  to  be  compared  to  the  need  to  solve  a  non-linear, 
non-convex  problem  in  the  classical  source  localization  prob¬ 
lem.  We  provide  a  proof  of  almost  sure  convergence  of  the 
MAP  solution  asymptotically  in  the  limit  of  a  large  num¬ 
ber  of  measurements,  a  step  that  tends  to  be  missing  in 
all  of  the  work  presented  earlier  in  the  source  localization 
literature. 

Second,  inspired  by  the  proof  of  convergence  of  the 
MAP  solution,  we  propose  and  implement  a  distributed 
sequential  regional  localization  algorithm:  Sense,  Trans¬ 
mit  &  Test.  This  algorithm  allows  for  sequential  sensing, 
transmission  and  testing  at  each  processor.  We  allow  each 
processor  to  have  one  or  multiple  regions  of  responsibility 
and  relate  the  probability  of  error  for  each  processor  in  the 
case  of  multiple  regions  to  the  probability  of  error  in  the 
case  of  a  single  region.  We  also  show  that  the  test  ends  in 
a  finite  time  under  mild  conditions  on  the  sensor  locations. 

Third,  we  illustrate  the  results  of  the  Sense ,  Transmit 
&  Test  and  show  how  the  expected  decision  time  for  a 
network  increases  with  the  required  accuracy  and  noise. 
We  also  provide  numerical  results  illustrating  how  it  is 
possible  to  increase  the  level  of  localization  accuracy  at 
the  expense  of  the  expected  decision  time  for  the  network 
for  a  fixed  decision  accuracy. 

l.f.  Paper  organization 

The  paper  proceeds  as  follows:  we  formulate  the  prob¬ 
lem  as  a  multi- hypotheses  testing  problem  in  Section  2. 
We  present  a  distributed  algorithm  to  solve  the  problem 
in  Section  3.  We  present  in  Section  4  numerical  results 
showing  the  performance  of  the  algorithm  as  various  pa¬ 
rameters  are  changed.  We  conclude  in  Section  5. 


distance.  The  nominal  distance  do  is  chosen  so  that  the 
received  power  in  the  vicinity  of  the  source  is  almost  equal 
to  the  transmitted  power  P  at  the  source.  Note  that  while 
this  model  gets  rid  of  the  singularity  at  the  source,  it  con¬ 
verges  to  the  same  behavior  as  the  classical  model  used  in 
communication  literature  Pi  =  \\q  Ps\\p^  when  the  distance 
|| qt  —  s||  is  large.  Here  P  is  the  power  received  at  a  unit 
distance  from  the  source.  The  received  power  becomes 

In  Pi  =  ln(Pdo)  -  ln(d0  +  || g*  -  s||p)  +  rq,  (1) 

where  n*  is  the  noise  associated  with  sensor  i,  and  all  ni 
are  independent  and  identically  distributed  (i.i.d)  Gaus¬ 
sian  random  variables  with  zero  mean  and  known  variance 
a2.  The  joint  probability  density  function  of  the  received 
power  Pr  =  [Pi,. . .  ,Pn]T ,  conditioned  on  the  source  lo¬ 
cation  y  €  Q  is 


p(P1,...,PN\y)  (2ircr2)N/2 


yN  , 

2 

A 

V  2a2  J 

(2) 


Problem  2.1  (MAP  point  localization  problem)  Com¬ 
pute  the  position  that  maximizes  the  conditional  density  of 
the  joint  observations,  that  is  compute 

y*  =  argmaxp(Pi, . . . ,  PN\y)V(y). 

2/6  <3 

Problem  2.1  is  a  nonlinear  nonconvex  optimization  prob¬ 
lem.  Attempts  to  solve  this  problem,  usually  revert  to 
relaxing  the  problem  or  approximating  its  solution  with¬ 
out  providing  a  convergence  analysis.  In  this  paper  we 
look  for  a  regional  localization,  so  the  conditioning  on  the 
exact  position  y  in  (2)  is  replaced  by  a  conditioning  on  the 
source  being  in  a  region  W*.  The  environment  Q  with  area 
A  is  divided  into  M  regions  { H'j , . . . ,  Wm}  with  positive 
areas  {Ai, . . . ,  Am}-  The  hypothesis  Ha  is  true  if  and  only 
if  s  €  Wa. 


2.  Source  localization  as  multi-hypothesis  testing 

We  start  this  section  by  introducing  the  model  and  the 
problem  definition. 

2.1.  Model  and  problem  definition 

Consider  a  compact  connected  environment  Q  C  R2. 
Suppose  that  there  are  N  sensors  placed  at  positions  qt  € 
Q  with  i  €  {1, . . . ,  A},  and  that  the  source  located  at 
an  unknown  location  s  £  Q  transmits  a  signal  whose 
power  undergoes  lognormal  shadowing  summarized  as  fol¬ 
lows.  The  average  power  loss  for  an  arbitrary  Transmitter- 
Receiver  separation  is  expressed  as  a  function  of  distance 
by  using  a  path  loss  exponent  p  >  2.  For  reasons  to  be 
explained  shortly,  we  work  with  a  slight  modification  of 
the  traditionally  used  model.  The  adopted  model  for  the 
received  power  at  a  sensor  i  is  Pi  =  do + a 1 1  p  >  where  p 
indicates  the  rate  at  which  the  power  loss  increases  with 


Problem  2.2  (MAP  regional  localization  problem) 

Compute  the  hypothesis  Ha  that  maximizes  the  posterior 
of  the  joint  observations,  that  is,  compute 

a*  =  argmax  p(P1, . . . ,  PN\Ha)¥(Ha).  (3) 

ae{l,...,M} 

2.2.  Regional  posterior  density 

Assuming  no  prior  knowledge  about  the  location  of  the 
source,  the  density  describing  s  €  Q  is 

=  /  l/A ’  if  s  G  Q’ 

^  I  0,  otherwise. 


Definition  2.3  (Repeated  measurements)  The  ith  sen¬ 
sor  takes  k  repeated  i.i.d.  noisy  measurements  and  com¬ 
putes  the  average  of  the  logarithms  of  the  measurements 


InP  i(k) 


E 


In  P»(Z) 
k 


(4) 


2 


In  the  infinite  measurement  case,  we  write 


In  P; 


In  Pi(t) 
k 


and  the  variance  lim^^  a2(k)  =  0. 

Proposition  2.4  (Expressions  for  posteriors)  In  the 

case  of  k  repeated  measurements,  the  regional  posterior  for 
sensor  i  about  region  Wa  is 


p(Pi(k)lHa)P(Ha) 


1/  _ I _ 

A  JWa  (27T(T2(fc))1/2 


x  exp  I 


(l"P<W-ln(A+fe||.) 

2  a2(k) 


dy, 


and  the  joint  regional  posterior  for  sensors  {1, . . . ,  N}  about 
region  Wa  is 


2.3.  Asymptotic  properties  of  regional  source  localization 
We  show  here  some  properties  of  the  MAP  algorithm 
when  applied  to  regional  source  localization  for  a  gen¬ 
eral  number  of  sensors  and  regions.  We  start  by  pre¬ 
senting  a  property  of  non-collinear  sensors  when  applied 
to  source  localization  using  measurements  undergoing  log¬ 
normal  shadowing. 


Lemma  2.5  (Three  non-collinear  sensors)  For  do  > 

0  and  p  >  0,  given  a  source  s  €  M2  and  three  non-collinear 
sensors  q±,  q2  and  q 3  €  M2,  the  only  solution  for  the  equa¬ 


tion 


EL  (: 


In 


4  I  .U  'a  |" 
<k>+i|s-g;|lp 


=  0  is  z  =  s. 


Proof:  In  fact,  it  is  easy  to  check  that  the  sum  is 
zero  at  z  =  s.  Uniqueness  of  this  solution  is  verified  by 
noting  that  the  sum  of  the  square  terms  is  zero  only  if  all 
the  summands  are  zero.  Let  q  =  (x,y)  and  qi  =  ( 5,1,  <7*2 )• 
The  solution  z  =  s  is  unique  if  and  only  if  the  following 
system  has  a  unique  solution: 


1  r 

-2(gn  -  <?2i) 

—2(512  —  522) 

X 

h 

p(Pi(k),...,PN(k)\HQ)P(Ha)  =  —  /  dy 

A  Jwa 

-2(511  -  531) 

-2(512  —  532)_ 

y_ 

k2 

A  1  ,  (hPiW-Wsrfc p)) 

(27ra2(fc))1/2eXPl  2cr2(fc) 

Proof:  Call  2  =  I11  P,(fc).  We  compute 


where  k\  and  k2  are  known  values  determined  by  the  mea¬ 
surements  and  the  positions  of  the  sensors.  The  system 
presented  in  Equation  (5)  has  a  unique  solution  if  and 
only  if  the  system  is  consistent  and  the  determinant  of  the 
matrix  is  non  zero,  i.e.,  the  three  points  are  non-collinear. 


,  ,  „  v  dP(Z<z,Ha)  d  fin  fWa  p(z\y)p(y)dydz 

P(2|HJ  =  Tz  r  W  =  Adz - K - 

_  d  flgp  fWa  (v(z\y) / A)dydz  _  fWa  p(zjy)dy 
dz  Aa  Aa 

Since  2  =  lnPj(fc)  =  J2t=i  the  conditional  proba¬ 

bility  is 

=  (27T(T2(fc))1/2 


As  usual,  assume  that  N  sensors  are  at  positions  qt,  i  £ 
{1,  ...,iV}  and  that  the  environment  is  partitioned  into 
closed  regions.  For  a  region  Wa ,  define  the  two  scalar 
quantities 


Ua 


L 


a. 


,  do+Wy-q.Y 

max  In  - - n - —  , 

yeWo,  d0+  s-5*  c 
JV} 


•  f-\  do  +  II V  ~  di llp 

mm  )  In  - - n - — 

v^Wa~i\  ^0  +  ||«  - '7i||p 


2 


(6) 

(7) 


x  exp  I 


(lnP.W-ln(jj+g»||,)) 


2a2 (k) 


dy. 


The  regional  posterior  is 


p(Pi(fc)|ffa)P(fl-a)  = 


Wa  (27T<T2(fc))1/2 

An 


x  exp  I 


2er2(fc) 


)dlJX 


fwa  (27r<T2(fe))1/2  'exp( 


2cr2(k) 


dy 


A 

Equations  for  the  joint  regional  posterior  follow  by  inde¬ 
pendence  of  measurements.  ■ 


Both  quantities  are  well  posed  because  they  are  the  max¬ 
imum  and  minimum  value  of  a  continuous  function  over  a 
compact  domain.  Additionally,  Ua  is  strictly  positive  for 
all  source  locations  s  £  Q  and  La  is  strictly  positive  for  all 
source  locations  s  £  Q  \  Wa .  The  latter  statement  follows 
from  Lemma  2.5  and  from  the  fact  that  the  distance  from 
s  to  Wa  is  strictly  positive  for  all  s  fL  Wa.  Define 

ya  =  \Ju2  +  ^-ua>  0,  (8) 

for  all  s  ^  Wa.  We  state  the  following  result  on  the  mag¬ 
nitude  of  sums  of  powers. 

Lemma  2.6  (On  the  posterior  of  a  wrong  hypothe¬ 
sis)  Consider  La,  Ua  and  ya  as  defined  in  (6),  (7)  and  (8). 
Assume  the  source  s  is  outside  Wa  and  the  noise  Ui  sat¬ 
isfies  \ni\  <  r]a  for  all  i  £  {1, . . . ,  N}  and  a  £  {1, . . . ,  M}. 
The  following  statements  hold: 


3 


1.  the  joint  measurement  is  lower  hounded  as 


Proof:  From  Lemma  2.6,  we  compute 


N 

min  y 

y£Wa  ^ 

X=1 


A  d0  +  \\y~qi\\P 
\  do  +  Ws-QiWP 


-  2L“’ 


and 


2.  the  posterior  probability  for  the  wrong  hypothesis  a 
is  upper  bounded  as 


p(P1,...,PN\Ha)F(Ha)  < 


Aa  exp(-La/4a2) 
A(2na2)N/2 


Proof:  To  prove  the  first  statement,  consider  the 

expansion 


d0  +  \\y-gi\\p 
V  do  +  \\s-qi\\P 


V  do  +  \\s  -  qi\\p  J 


+  2 


do  +  II y  - 

do  +  \\s  -  qi\\p  )  1  f-f 


By  computing  lower  bounds  for  each  term  and  substituting 
the  definition  of  ya,  obtain 


N 

min  y  (  In 

y&Wa  \ 


\\y-Qi\ 


,  ,  I,  \\p+ni)  ^  La-2NUar]a-Nr]l 
d0  +  \\s-qi\\p  ) 

=  La+  2NU2a  -  2NUa]ju2  +  ^-  N(ul  + 


-NUl  +  2NUaXIUl  +  ±A=l-L 


P[p(P1,...,PN\Ha)V(Ha)<ea(a)\ 
>P[[n1,...,nJV]T  G  [~Va,Va]N] 

N  /f  1  \ 

=  n  (^2  P[ni  >  da\  +  2  “  P[n*  <  ~rl‘x}) 

=  (l  -  2Q(?7a/o-))Ar. 

The  first  inequality  follows  from  the  fact  that  Lemma  2.6 
holds  whenever  all  n,;  <  r]a.  The  proofs  of  the  two  limits 
of  limfc^oo  ea(tXfc)  and  lim^oo  lia(ak)  are  immediate  when 
there  are  at  least  3  non-collinear  sensors.  Indeed,  if  there 
are  at  least  3  non-collinear  sensors  and  if  s  Wa,  then 
Lemma  2.5  applies  and  one  can  show  La  >  0  and  rja  >  0. 

■ 

This  theorem  states  that,  as  a  — >  0+,  the  joint  regional 
posterior  p(Pl5 . . . ,  P/v|iLQ)P(iLQ)  takes  an  arbitrarily  small 
value  with  a  probability  that  goes  arbitrarily  close  to  1 
when  Ha  is  not  the  correct  hypothesis.  This  is  so  as 
Q(x )  — »  0  as  x  — >  oo.  To  complement  the  Theorem  2.7, 
we  prove  below  that  for  the  correct  hypothesis,  the  prob¬ 
ability  density  is  lower  bounded  by  a  positive  term  w.p.l. 

Theorem  2.8  (Strict  positivity  of  correct  hypothe¬ 
sis)  Consider  sensors  at  positions  qi, . . . ,  qjy.  Let  a  be  the 
noise  variance.  If  the  source  s  G  Wg,  then 


P  [p(P1;  .  .  .  ,  PN\Ha)P(Ha)  >  ’F(ct)]  >  n(<7), 


The  second  statement  follows  directly  from  the  first  state¬ 
ment  because  of  the  equality 


In  Pt  —  In 


Pd0 

d0  +  II y  -  <ti\\p 


,  do  +  ||y  -  qi\\p 
do  +  \\s-qi\\P 


+  rii, 


and  because  of  the  fact  that  the  surface  integral  of  a  func¬ 
tion  /  is  upper  bounded  by  the  surface  integral  of  the 
maximum  value  of  /.  ■ 

We  are  now  ready  for  the  convergence  theorem.  We  intro¬ 
duce  the  standard  function  Q  :  K.  — >  K>0  by 

1  f+°° 

Q{%)  =  /  exp(~y2/2)dy. 

V  "TT  J  X 


Theorem  2.7  (Elimination  of  wrong  hypothesis)  Con¬ 
sider  sensors  at  positions  q±, . . . ,  <pv-  Let  a  be  the  noise 
variance.  If  the  source  s  £  Wa,  then 


p(Pi, . . .  ,PN\Ha)P(Ha)  <  ea(cr) 


>  Ha(cr), 


where 

Aaexp(-La/4a2)  N 

tocW)  = — A(flTTa2)N/2 — ’  “  2(3(??a/cr))  • 


where 


T(cr)  =  p(Pi,...,Pjv) 
fl(cr)  =  /J.a(cr) 

at^a 


-  E 

a=l 

OL^Ot 


Aa  exp(— La/4tJ2) 
A(2na2)N I2 


=  II  (1-2  Q(Va/a))N. 

oc=l 


Furthermore,  in  the  k  repeated  measurement  case,  if  at 
least  3  sensors  are  non-collinear,  then  lim^oo  (crfc)  = 
p(Pi ,  •  •  • ,  Pn)  >  0  and  lim^oo  fl(ak)  =  1~ . 

Proof:  The  proof  of  this  theorem  follows  directly 

from  Theorem  2.7  and  from  the  total  probability  theorem. 
Call  0  =  [Pi, ... ,  Pn)t.  We  know  from  the  total  probabil¬ 
ity  theorem  that 

M 

P(*0  =  ^2p{z\Ha)P(Ha)  =  p(z\Hg)P(Hg) 

Ct=  1 

+  Y  P(z\Ha)V(Ha) 


and,  in  turn,  that 


Furthermore,  in  the  k  repeated  measurement  case,  if  at 
least  3  sensors  are  non-collinear,  then  lim^oo  ea(ak)  = 
0+  and  lim^oo  pa(ak)  —  1~ . 


p(z\Hw)¥(Hw)  =  p(2)  -  Y  P(2I Ha)P(Ha). 

Oiz^a. 
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From  Theorem  2.7 


p(z\Ha)V(Ha)  >  PO) 


ac=l,...,M 

OL^a. 


>  n 


p{z\Ha)P(Ha)  <  ea(a) 


>  n  M0") 


As  limfc^oo  CTfc  =  0+,  limfc^oo  'F(crfc)  =  p(z)  and 
limfc_»oo  ^(cfc)  =  l--  ■ 

Theorem  2.8  complements  Theorem  2.7  in  that  is  shows 
that  as  lim/.^oo  <jk  =  0+,  the  largest  regional  posterior  is 
the  one  associated  with  the  correct  hypothesis. 


Remark  2.9  (Almost  sure  convergence  of  MAP)  Us¬ 
ing  a  MAP  algorithm  to  solve  the  problem  of  regional  lo¬ 
calization,  is  assured  to  provide  a  correct  answer,  almost 
surely,  in  the  limit  of  infinite  measurements.  This  follows 
directly  from  Theorem  2.8  and  Theorem  2.7. 


3.  Distributed  sequential  regional  localization 

In  this  section  we  assume  that  each  sensor  is  a  processor 
that  can  perform  computational  tasks  as  well  as  commu¬ 
nicate  to  other  processors  according  to  a  specified  commu¬ 
nication  graph.  Each  processor  takes  measurements  and 
computes  a  conditional  posterior  that  it  communicates  to 
all  its  neighbors  and  then  makes  a  decision  if  a  desired 
accuracy  is  reached.  A  group  of  regions  is  associated  with 
each  processor.  The  processor  will  need  to  provide  a  de¬ 
cision  about  which  of  these  regions  if  any  contains  the 
source.  We  call  such  a  group,  the  regions  of  responsibility 
of  the  processor.  We  do  not  assume  any  constraints  on  the 
assignment  of  regions  of  responsibilities.  We  present  the 
algorithm  in  Subsection  3.1  and  describe  its  properties  in 
Subsection  3.2. 

3.1.  Distributed  algorithm  based  on  sequential  sensing,  com¬ 
munication  and  hypothesis  testing 

We  present  below  a  distributed  algorithm  where  each 
processor  decides  whether  or  not  its  region  of  responsi¬ 
bility  contains  the  source.  The  algorithm  as  presented, 
dictates  until  when  a  processor  needs  to  continue  to  take 
measurements,  as  well  as  the  information  that  needs  to  be 
communicated. 

For  each  processor  i  £  {1, . . . ,  TV},  the  set  of  neighbors 
Ni  consists  of  the  processor  itself  along  with  the  proces¬ 
sors  that  can  communicate  with  it.  The  it li  processor  is 
responsible  for  a  set  lZt  of  Mi  regions.  We  denote  these 
Mi  regions  by  Wa  for  a  £  TZ  {.  The  processor  collects  the 
measurements  from  its  neighboring  processors,  and  calcu¬ 
lates  two  posteriors  for  all  regions  Wa,  a  £  IZi-  The  first 
posterior  corresponds  to  the  hypothesis  that  the  source  is 
in  W a ,  the  second  posterior  corresponds  to  the  hypoth¬ 
esis  that  the  source  is  outside  Wa.  Once  the  processor 


reaches  a  pre-defined  level  of  confidence,  it  provides  a  de¬ 
cision  about  whether  or  not  Wa  contains  the  source.  The 
itli  processor  stops  running  its  test  when  it  reaches  a  de¬ 
cision  about  all  Wa,  a  £  IZi.  The  processor  then  sets  its 
decision  to  either  yes,  the  source  is  in  Wa,  or  no,  no  source 
is  in  Gaem  Wa.  Each  processor  continues  to  sense  and 
transmit  its  measurements  until  all  its  neighbors  j  £  Mi 
have  reached  a  decision.  We  give  here  a  formal  description 
of  the  algorithm. 

Algorithm  :  Sense,  Transmit  &  Test 

algorithm  tolerance:  0  <  e<  ) 
network  processors:  *  £  {1, . . . ,  N} 
regions:  WA,  a  £  {1 , ...  ,M} 
state  of  processor  i  contains: 

a-dcsn;  £  {yes  source  £  WA,no  source  £  LA  WA, unknown}, 
for  all  j  £  Mi  :  Q.;,a-stop(  £  {false,  true}, 
for  all  a  £  IZi  '■  WA,r-stopa  £  {false,  true}, 
r-dcsnQ  £  {yes, no,  unknown} 

Processor  i  with  set  of  neighbors  Mi  executes: 

1:  transmit  qi  to  j  £  Mi 

2:  set  k  :=  0  ,  a-stop,  :=  false  and  a-dcsn;  :=  unknown 
3:  set  r-stopa  :=  false  and  r-dcsnQ  :=  unknown  for  a  £  IZi 
4:  While  3  j  £  Mi  with  a-stop^  ==  false  do 
5:  update  k  :=  k  +  1  and  take  measurement  Pi(k) 

6:  compute  lnP;(fc)  =  X)t=i  lnt^^ 

7:  transmit  lnP;(fc)  to  j  £  Mi 

8:  store  lnPjv^fc)  =  {lnP;(fc)}  U{lnPj(fe)  for  all  j  £  Mi} 

9:  For  all  a  £  IZi  with  r-stopa  — =  false  do 

10:  If  p(Pjwi(*)|s  €  WA)P(s  €  WA)  >  (1  -  e)  p(Pjv4(fc)) 

11:  dcsnQ  :=  yes,  r-stopa  “true,  a-dcsn;  :=  true 

12:  If  p(P^(k)\s  i  WA)P(£  WA)  >  (1  -  e)  p(PAf,(fc)) 

13:  r-dcsn^  :=  no  and  r-stopa  :=  true 

14:  End  For 

15:  If  r-stopQ  ==  true  for  all  a  £  IZi 

16:  a-stopi  :=  true 

17:  If  a-dcsn;  ==  unknown 

18:  a-dcsn,  :=  no 

19:  transmit  a-stop;  to  all  j  £  Mi 

20:  return  a-dcsn, 

21:  End  While 


3.2.  Properties  of  Sense,  Transmit  &  Test 

We  present  below  properties  involving  the  accuracy  and 
decision  time  of  the  Sense,  Transmit  &  Test  algorithm. 

Theorem  3.1  (Accuracy  and  decision  time  for  Sense, 
Transmit  &  Test  algorithm)  Assume  that  only  one  source 
exists  in  the  environment  Q,  that  each  processor  has  at  least 
2  neighboring  processors  with  which  it  forms  a  non-collinear 
triplet,  and  that  each  processor  is  assigned  Mi  regions.  Given 
an  accuracy  e  £  (0,  |),  the  Sense,  Transmit  &  Test  algorithm 
enjoys  the  following  two  properties: 

1.  the  algorithm  ends  in  a  finite  time,  and 

2.  each  processor  i  has  a  probability  of  error  no  larger  than 
2 Mit  if  2  <  Mi  <  1  +  i,  and  no  larger  than  e  if  Mi  =  1. 

Proof:  It  is  well  known  (Wald,  1945;  Varshney,  1996)  that 
given  two  hypothesis  Hi  and  Ho  with  known  posteriors,  P(Hi) 
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and  P(H0),  a  hypothesis  test  that  ensures  that  the  decision  un¬ 
der  hypotheses  Ho  and  Hi  is  correct  with  a  probability  greater 
than  to  and  n  respectively  is  the  following: 


1. 

2. 

3. 

4. 


Calculate  p(PArj(A:)|fl'i)P(//1 

p(PAf.(fc)|,F/l)P(ff1) 


if 


if 


_  >  T1 

Pip^.(t)iffo)nHo)  -  i-Ti 
piPA/vWiffimwi)  ^  i-T-o 
1-0 


,  p(PjV;  {k)\Ho)V(Ho) 
decide  in  favor  of  Hi, 


< 


decide  in  favor  of  Hi, 


p(P^(fc)|ffo)P(»o)  - 
othewise  repeat  measurements  and  go  to  1. 


We  show  below  that  the  Sense,  Transmit  &  Test  algorithm 
satisfies  the  description  above.  Applying  the  total  probability 
theorem  we  get 


p{P#t\8$Wa)P(8$Wa) 

=  P (Put  | y  €  Q)P(s  €  Q)  -  p(PatJs  €  Wa)V(s  €  Wa) 

=  p(JV1)-p(iVi|a6W’a)P(s6W’„).  (9) 

Call  (Hi  :=s£  Wa)  and  (Ho  :=  s  ^  Wa).  If  we  set  to  =  n  = 
(1  —  t),  and  the  thresholds  to  accept  a  hypothesis  Hi,  to  be 


p(PMi(k)\HiMHi)  >  rip (PMi), 


and  the  thresholds  to  reject  a  hypothesis  to  be 


To  write  a  formal  proof,  we  introduce  Pf  and  pm  to  be  the 
probability  of  false  alarm  and  mis-detection.  Where  pf  cor¬ 
responds  to  choosing  yes  while  the  correct  decision  is  no  and 
mis-detection  corresponds  to  choose  no  when  the  correct  de¬ 
cision  is  yes  for  any  region  Wa.  Here  that  Pf  =  Pm  =  t.  A 
processor  makes  an  error  if  it  wrongly  decides  yes/no  on  Wa 
for  any  a  £  Hi.  Following  the  analysis  above,  the  probability 
of  error  for  the  ith  processor  is: 

Pe  <  [pmP(s  i  Ua6TC,  Wa) 

+  (Pf  +  ^jPfP^)P(S  £  Ua6R,  Wa))  <  2 Mit, 

if  e(Mi  —  1)  <  1.  If  the  processor  has  only  one  region  of  re¬ 
sponsibility,  it  is  straightforward  to  see  that  the  processor  has 
a  probability  of  error  no  larger  than  e. 

We  show  now  that  the  test  ends  after  a  finite  number  of 
measurements.  For  a  region  Wa,  processor  needs  to  decide 
whether  the  source  is  in  Wa  (Hi)  or  outside  it  (Ho). 

Without  loss  of  generality,  assume  that  Hi  is  correct  for  a 
region  W„, .  We  know  from  Theorem  2.7  that 


p(PuM\HomHo)  >  top(Pm), 
then  using  Eq.  (9),  one  can  show  that 

p(PJvi(fc)|P1)P(Jff1)  >  n  p (PMi) 
p(PAri(fe)|Po)P(Po)  <  (1  -  Ti)p(PMi) 
p(P,v,(fc)|Pi)P(Pi)  n 
(k)\H0)F(H0)  ~  1-Ti 

Similarly,  assuming  Ho  is  correct,  one  can  show  that 


lim  P(Pa7  (k)\Ho)¥(Ho)  =  0+, 

k — >oo 

almost  surely.  We  also  know  from  Theorem  2.8  that 
lim  p(PMi(k)\Hi)F(Hi)  =  p(Pv4(fc))  >  0, 

k — >oo 

almost  surely.  This  has  the  following  implication 

hm  P(p^wi^)n^)  =  lf 

fe^OO  p(P//((«)) 


p(PjVz(k)\Ho)¥(Ho)  >  to  p(PMi) 

p(PMl(k)\Hi)V(Hi)  <1~rp 
p(PMi(k)\Ho)¥(H0)  -to  '  V  j 

Assuming  Hi  is  correct,  the  probability  of  correct  decision  for 
the  ith  processor  is  no  smaller  than  n,  for  each  of  the  regions 
Wa,ot  £  Mi.  Similar  result  hold  assuming  Ho  is  correct. 

The  maximum  numbers  of  errors  that  a  processor  can  make 
in  a  decision  is  two:  a  mis-detection  and  a  false-alarm  for  Wa 
where  a  £  Hi.  Alternatively  all  the  other  combinations  of 
choices  result  in  at  most  one  error,  since  the  *th  processor  can 
declare  at  most  one  hypothesis  Ha  to  be  correct  for  all  a  £  Hi. 

The  scenarios  where  the  decision  of  the  processor  is  erro¬ 
neous  are  presented  below: 

1.  If  one  of  the  regions  Wa  satisfies  p(PJy’j  (k)\Hi)¥(Hi)  > 
Tip(PUi),  then  for  all  (3  £  Hi  \  a,  the  following  holds 
(from  the  complete  probability  theorem) 

p(PMi(k)\s  £  Wp)P(s  £  Wp)  <  (1  -  n)p(JV4) 

=»p(PAAi(fc)|s  ^  Wp)f(s  £  Wp)  >  Tl  P (PUi)  =  top(Pat,). 

(12) 

From  Eqs.  (11)  and  (12)  it  follows  that  the  source  can 
be  detected  in  at  most  one  region  Wa.  It  follows  that  at 
most  one  false  alarm  can  happen,  which  might  or  might 
not  be  accompanied  with  one  mis-detection. 

2.  If  none  of  the  regions  of  responsibilities  of  the  processor 
contain  the  source,  then  the  processor  can  make  at  most 
one  mistake  by  having  at  most  one  false  alarm. 


which  implies  that  for  all  e  >  0,  there  exists  0  <  K  <  oo,  s.t. 

,  P(PW (K)\Hi)F(Hi)  -  p(P^  (A'))  „  _ 

p(Pv,W)  1 

-  p(PJvl(A')|Hi)P(H1)  +  p(PArj(AT))  <  e  p(P^(I<)) 
p(PM,(K)\Hi)¥(Hi)  >  (1  -  e)  p(PK(K)). 

So  for  any  |  <  r  <  1,  there  exists,  almost  surely,  K  <  oo,  s.t. 

p(PUi(k)\Hi)V(Hi)  >  r  p(PM(fc)), 

here  r  =  1  —  e,  where  0  <  1  <  \ . 

Similarly  one  can  prove  that  if  Ho  is  correct,  then  there 
exists,  almost  surely,  K  <  oo,  such  that 

v(Puz(k)\Ho)^(Ho)  >  t  p(Pjs/i(fc)). 

To  complete  the  proof,  we  cover  the  cases  where  the  algo¬ 
rithm  makes  a  wrong  decision.  This  is  possible  if  the  thresh¬ 
olds  corresponding  to  a  wrong  decision  are  crossed  at  a  time 
K\  <  K  <  oo. 

This  completes  the  proof  that  the  Sense,  Transmit  &  Test 
algorithm  has  a  finite  decision  time.  ■ 

4.  Numerical  results 

We  present  in  this  section  three  sets  of  simulations.  The 
first  two  sets  illustrate  some  properties  of  the  Sense,  Transmit 
&  Test  algorithm,  while  the  third  presents  a  modification  of 
the  algorithm  that  introduces  an  interesting  extension  of  the 
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work.  In  the  first  simulations,  there  are  as  many  regions  as 
there  are  sensors,  i.e. ,  N  =  M  =  10.  We  start  by  presenting  in 
Figure  1  a  sample  of  the  results  obtained  by  the  Sense,  Trans¬ 
mit  &  Test  algorithm.  The  figure  shows  the  positions  of  the 
processors  (equipped  with  sensors)  as  well  as  the  partition  of 
Q.  As  partition  we  adopt  the  Voronoi  partition  generated  by 
the  processors  positions;  each  processor  is  responsible  for  its 
corresponding  Voronoi  region.  As  stated  in  the  caption,  after 
113  observations  all  decisions  have  been  made  and  the  source 
has  been  correctly  localized. 


Figure  1:  This  picture  illustrates  an  evolution  of  the  output  of  the 
Sense ,  Transmit  &  Test  algorithm.  At  each  instant  a  region  is  col¬ 
ored  in  white,  light  gray  or  dark  gray,  indicating  unknown,  yes  or  no 
respectively.  The  output  of  the  distributed  algorithm  is  shown  at 
times  0, 1,  3,  4,  6,  8, 11, 113  respectively.  In  this  run  we  e  =  0.01  and 
a  =  0.5  with  N  =  M  =  10. 


We  then  present  in  Figure  2  a  plot  that  shows  how  the  ex¬ 
pected  number  of  observations  needed  to  reach  a  decision  varies 
with  the  accuracy  t  in  the  algorithm.  Clearly,  the  probability 
of  correct  detection  increases  for  decreasing  values  of  e. 
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Figure  2:  This  plot  shows  the  expected  time  it  takes  a  network  of 
10  processors  implementing  the  Sense,  Transmit  &  Test  algorithm 
to  reach  a  decision  for  a  noise  standard  of  deviation  a  =  0.1  when 
the  probability  of  error  e  varies.  We  show  the  logarithm  of  the  de¬ 
cision  time.  Note  that  the  network  decision  time  seems  to  grow 
exponentially  with  the  desired  accuracy  as  is  standard  in  sequential 
hypothesis  testing.  The  network  is  assumed  to  have  reached  a  deci¬ 
sion  when  all  processors  have  decided.  The  expected  decision  time 
is  calculated  over  1000  runs. 

In  Figure  3  we  show  how  the  expected  number  of  obser¬ 
vations  needed  to  reach  a  decision  increases  with  the  standard 
deviation  of  the  noise. 

Next,  we  report  the  second  sets  of  simulations,  where  we 
have  differing  numbers  of  regions  and  sensors.  Specifically,  we 
have  N  =  4  sensors  and  M  =  16  regions.  Figure  4  illustrates 
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Figure  3:  This  plot  shows  the  expected  time  it  takes  a  network  of  10 
processors  implementing  the  Sense,  Transmit  &  Test  algorithm  to 
reach  a  decision  with  a  probability  of  error  no  larger  than  e  =  0.01 
as  the  noise  standard  of  deviation  <r  varies.  The  network  is  assumed 
to  have  reached  a  decision  when  all  processors  have  decided.  The 
expected  decision  time  is  calculated  over  100  runs. 


the  evolution  of  the  Sense,  Transmit  &  Test  algorithm  in  this 
case.  The  overall  accuracy  for  each  processor  is  0.9.  This  is 
achieved  by  setting  e  =  0.1/8. 


Figure  4:  This  picture  illustrates  an  evolution  of  the  Sense,  Trans¬ 
mit  &  Test  algorithm.  At  each  instant  a  region  is  colored  in 
white,  light  gray  or  dark  gray,  indicating  unknown,  yes  or  no  respec¬ 
tively.  The  output  of  the  distributed  algorithm  is  shown  at  times 
1,4,6,7,13,131,142,202  respectively.  In  this  run  we  set  e  =  0.1/8 
and  a  =  0.5  and  N  =  4  and  M  =  16. 

In  this  third  set  of  simulations,  we  show  the  output  of  a 
modified,  multi-resolution  version  of  the  Sense,  Transmit  & 
Test  algorithm.  This  multi-resolution  version  runs  over  multi¬ 
ple  stages,  at  each  stage  the  environment  under  consideration 
is  divided  in  two  regions.  Observations  are  taken  at  each  stage 
until  one  of  the  two  regions  is  rejected  with  an  accuracy  of  1  —  e. 
The  rejected  region  is  removed  from  the  environment,  and  the 
remaining  region  is  again  divided  in  two  regions.  Observations 
are  transferred  from  one  stage  to  another  and  re-used  to  reach 
a  decision  about  the  more  fine  environment  division.  A  sample 
output  of  the  modified  algorithm  is  shown  in  Figure  5.  In  or¬ 
der  to  reach  the  same  precision  in  localization  as  that  shown  in 
Figure  4,  we  divide  the  regions  4  times.  Note  that  the  original 
Sense,  Transmit  &  Test  algorithm  reached  its  decision  after 
an  average  of  290  observations  whereas  the  multi-resolution 
algorithm  did  so  after  an  average  of  100  observations.  We  cal¬ 
culated  these  values  from  1000  Monte-Carlo  runs,  that  is  with 
an  error  of  ±3%  to  show  a  similar  probability  of  error  with 
the  same  level  of  fineness.  We  leave  a  rigorous  analysis  of  the 
multi-resolution  algorithm  to  future  work. 
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Figure  5:  This  picture  illustrates  an  evolution  of  the  modified  version 
of  the  Sense,  Transmit  &  Test  algorithm.  The  output  is  shown  at 
times  7,11,71,72  respectively.  The  rejected  regions  are  shown  in 
dark  grey  and  the  ones  accepted  at  each  stage  are  shown  in  light 
grey.  In  this  run  we  set  e  =  0.1  and  <r  =  0.5  and  N  =  4  and  M  =  2 
at  each  set  of  tests. 

We  conclude  with  a  general  remark.  The  Sense,  Transmit  & 
Test  algorithm  presented  in  this  work  might  at  first  glance  look 
similar  to  sequential  multiple  hypothesis  testing  algorithms  by 
elimination  such  as  the  one  presented  in  (Bauer,  1989).  A  closer 
comparison  of  the  two  algorithms  shows  that  while  in  this  work 
at  most  2 M  tests  are  run  at  each  sample,  the  hypothesis  test  by 
sequential  elimination  requires  a  number  of  tests  of  the  order 
2m  as  it  proceeds  by  a  pairwise  comparison  over  all  hypothesis. 
Nonetheless,  it  is  worth  mentioning  that  while  the  sequential 
elimination  algorithm  leads  to  a  decision  as  soon  as  all  but  one 
hypothesis  is  eliminated,  we  wait  here  until  the  last  hypothesis 
reaches  the  required  certainty  level.  This  can  be  seen  in  Fig¬ 
ure  1  where  all  but  one  hypothesis  were  canceled  at  the  11th 
observation,  yet  the  algorithm  did  not  end  until  the  113th  ob¬ 
servation  when  the  last  processor  reached  its  required  accuracy. 
The  geometric  aspects  and  the  properties  associated  with  the 
regional  localization  problem,  made  it  possible  to  propose  the 
simpler,  yet  less  general,  Sense,  Transmit  &  Test  algorithm. 

5.  Conclusion 

In  this  work,  we  looked  at  the  problem  of  source  localiza¬ 
tion  in  a  multiple  hypothesis  testing  setting.  We  based  our 
formulation  on  the  geometric  properties  of  the  MAP  algorithm 
when  applied  to  regional  localization.  We  proved  that  when 
measurements  are  available  from  three  or  more  non-collinear 
sensors,  MAP  based  algorithms  choose  the  correct  region  al¬ 
most  surely  in  the  limit  of  infinite  measurements.  We  then 
presented  a  sequential  distributed  algorithm  where  each  pro¬ 
cessor  senses,  transmits  and  tests  to  provide  a  decision.  We 
analyzed  the  algorithm  and  provided  a  measure  of  its  accuracy 
and  showed  that  it  ends  in  a  finite  time.  We  concluded  the 
paper  by  numerically  illustrating  the  algorithm’s  performance. 

There  are  two  direct  extensions  for  this  work  that  we  are 
considering.  The  first  is  using  an  adaptive  hierarchical  methods 
based  on  quadtrees  (dc  Berg  et  al.,  2000)  to  increase  the  level 
of  details  in  the  choice  of  regions.  The  regions  could  be  finely 
divided  as  fewer  candidate  regions  are  left,  an  example  of  such 
adaptation  is  shown  at  the  end  of  the  manuscript.  It  would  be 
interesting  to  study  the  trade  off  between  the  accuracy  and  the 
decision  time  as  a  function  of  the  fine-gridding  of  the  regions. 

The  second  is  allowing  the  algorithm  to  stop  as  soon  as  a 
processor  decides  that  its  region  contains  the  source.  As  pre¬ 
sented  in  this  manuscript,  the  algorithm  has  a  proven  accuracy 
performance  based  on  the  assumption  that  all  processors  reach 
their  decisions  independently  of  each  other,  and  although  we 
assume  only  one  source,  a  processor  will  continue  applying  the 
Sense,  Transmit  &  Test  algorithm  until  it  decides  that  its  re¬ 
gion  does  not  contain  the  source  even  if  the  source  was  detected 
by  one  of  the  other  processors.  It  will  be  interesting  to  see  what 


happens  to  the  accuracy  if  an  individual  can  broadcast  a  yes  to 
everyone  in  the  group,  allowing  them  to  stop.  Alternatively,  as 
we  showed  in  Figure  1,  it  is  possible  that  only  one  hypothesis 
is  left  by  elimination.  It  will  also  be  interesting  to  analyze,  if 
possible,  the  accuracy  of  an  algorithm  that  makes  use  of  such 
scenarios  when  they  occur. 
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